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Patterns of vortex ripples form when a sand bed is subjected to an oscillatory fluid flow. Here 
we describe experiments on the response of regular vortex ripple patterns to sudden changes of the 
driving amplitude a or frequency /. A sufficient decrease of / leads to a "freezing" of the pattern, 
while a sufficient increase of / leads to a supercritical secondary "pearling" instability. Sufflcient 
changes in the amplitude a lead to subcritical secondary "doubling" and "bulging" instabilities. Our 
findings are summarized in a "stability balloon" for vortex ripple pattern formation. 
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A flat bed of sand subjected to an oscillatory flow of 
water is seldom stable but instead displays the forma- 
tion of patterns. Classical studies (l) have shown that 
after the flat bed becomes linearly unstable, so-called 
rolling grain ripples (small heaps of grains) are formed 
first. These are, however, always transient l2l, and even- 
tually strongly nonlinear vortex ripples are formed in a 
coarsening type process. These ripples have triangular 
crests with slopes roughly at the angle of repose, and the 
flow around the ripple crests is dominated by vortices 
that occur in the wake of the ripples. Together with the 
converging flow at the "wind" side of the ripples, these 
vortices yield sand mass transport directed toward the 
crests of the ripples, which is balanced by sand avalanch- 
ing down when the slopes grow too large. The wavelength 
of such ripple patterns is comparable to the amplitude of 
the fluid motion which sets the scale for the size of the 
separation bubbles ; this wavelength is substantially 
larger than the most unstable wavelength of a flat bed 

i|. 

To characterize the vortex ripples we have studied their 
pattern forming properties. From this perspective, the 
system combines a number of unique features. Firstly, 
the driving is anisotropic which results in alignment of 
the ripples perpendicular to the flow. This allows for 
studies in one-dimensional geometries |^J^ , although our 
experiments indicate that instabilities of the ripples lead 
to intrinsically two dimensional patterns. Secondly, typ- 
ical ripple wavelengths are essentially independent of the 
system and grain dimensions and /, but scale with the 
driving amplitude a Finally, due to the strongly 

nonlinear character of the development, it has so far not 
been possible to describe the pattern dynamics in terms 
of an "amplitude equation" . 

Our setup consists of a tray of sand that is oscillated 
with amplitude a and frequency / in a tank of water, 
allowing us to study two-dimensional patterns. To probe 
the fully developed vortex ripples, we have studied their 
response to sudden changes of the control parameters a 
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FIG. 1. Stability balloon for fully developed vortex rip- 
ples. The initial ripple pattern with wavelength 4.2 cm is 
compatible with driving parameters a = 2.8 cm and / = 0.6 Hz 
(cross) and remains stable for drivings indicated in full circles. 
Triangles, squares and diamonds indicate drivings where this 
pattern experiences a doubling, pearling or bulging instabil- 
ity. The open circles correspond to the n = 10 measurements 
shown in Fig. ^ and the dashed line roughly indicates where 
the pattern freezes. 

and /. Our flndings lead to a "stability balloon" for 
vortex ripple pattern formation shown in Fig. |l|. First 
we study the amount of grain motion for patterns that 
have evolved "freely" from the flat bed and have a wave- 
length selected by the driving amplitude. This leads to 
the "freezing line" in Fig. |l| (and discussed in more detail 
in Fig. ^). In the second set of experiments we study 
the response of regular patterns with initially fixed wave- 
lengths to changes in / or a. We find that a secondary 
"pearling" instability (Fig. ^ occurs when the frequency 
is sufficiently increased, while secondary "doubling" 
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FIG. 2. Sketch of the experimental setup. The tray A is 
filled with sand and its oscillations are driven by the motor F 
via the rail B. 

and "bulging" (Fig. ^ instabilities occur when the am- 
plitude a is decreased or increased respectively, beyond 
some threshold. A qualitative presentation of these in- 
stabilities has appeared in 

Experimental Setup ~ Our setup is sketched in Fig. 
An 0.6 mx 1 m aluminum tray (A) suspended from rollers 
(B) mounted on a stiff frame (C) is immersed in a suffi- 
ciently deep tank (D) filled with water. To avoid slosh- 
ing, a 3 cm thick flat plate of Plexiglas (E) is placed on 
top of the tank. The frame and tray are driven with 
a continuously controllable frequency / (period T) and 
amphtude a by an ICME AC motor (F) . The sides of the 
tray consist of 2 cm high straight boundaries, while the 
ends are triangular wedges with a slope of 15° and max- 
imum height 2 cm (see Fig. ^ ; the rational behind these 
"soft" boundaries will be discussed below. The "sand" 
consists of spherical glass beads, ranging in sizes from 
250 to 350 /im. The thickness of the sand layer is small- 
est in the troughs of the ripple pattern, but always larger 
than 5 mm. The whole setup is illuminated from the left 
side and filmed from above by a Dalsa 8 bit CAD4 CCD 
camera with 1024^ resolution. A trigger is mounted on 
the motor so that all pictures are taken at the same ex- 
tremal position of the tray. For the values of the driving 
considered here (Fig. |^), typical ripple patterns consists 
of 15-20 ripple lengths. Suspension can be ignored and 
the maximal acceleration of the tray is well below the 
fluidization threshold. For appropriate values of a and 
/, fairly regular ripple patterns grow from the flat bed. 
For example, for a = 2.8 cm we find ripple patterns with 
wavelength 4.2 cm. 

Grain motion - Once a fully developed pattern is 
formed, how does the number of grains in motion vary 
with / and a? We observed that, due to irregularities in 
some grains, sand bed images usually display a number of 
very bright spots. The difference between two subsequent 
images is dominated by the appearance or disappearance 
of a number n of such bright spots; we assume that n is 
proportional to the number of grains which have moved. 

To measure n we proceed as follows: We start from a 
ffat bed and obtain an equilibrated pattern by running 
the system for one hour, so that the driving amplitude 
selects the ripple wavelength, in contrast to the secondary 
instability experiments presented below. We then take 
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FIG. 3. (a) The average amount of moving grains n as 
a function of the driving frequency / for amplitudes 1.4 
(dot-dashed), 2.4 (dashed) and 4.4 (dots), (b) Points in con- 
trol parameter space where n is 10 (open circles) and 50 
(closed circles) compared to a curve where aJ^'^p is con- 
stant, (c-d) Data collapse of ^Jn versus the turbulent ex- 
pression a'^^f^ and the laminar expression Ap^^ , with a in 
cm and / in Hz (see p^). 

a series 50 images from which the average of n is de- 
termined. The frequency is then lowered to a new value, 
the system is allowed to relax for fifteen minutes and fifty 
new images are recorded. This procedure is repeated for 
ten decreasing values of the frequency. In Fig. ||a, three 
examples of the number of moving grains as function of 
frequency are shown. For high frequencies n saturates, 
and the corresponding data points were discarded in the 
analysis presented below; this leaves us with a total of 65 
data points for 10 different values of a. 

We wish to obtain a relation between a and / where 
n is constant. It was previously suggested that the max- 
imum non-dimensional shear stress on a flat bed, 
would be a good measure of the forcing |,|,|l3- If 
the boundary layer were laminar one would obtain that 
a-m oc aj"^!^ from the solution to Stokes second problem 
pT[ , while for turbulent flow a semi-empirical relation 
am oc a^'^^/^ would be appropriate jl^. In Fig. 0c we 
have plotted \pn versus the turbulent expression a^'"^^ f"^ 
and obtain a fairly good linear relation, while the cor- 
relation between n and the laminar expression a]"^!^ is 
very weak (Fig. |d). 

Secondary instabilities - So far we have described pat- 
terns with a wavelength that is selected by the driving 
amplitude. Now we ask what happens when a perfectly 
regular pattern is suddenly subjected to changes in the 
driving parameters a or /. 

Obtaining a completely regular pattern is not entirely 
trivial. Freely grown patterns contain defects, which may 
be annihilated only after long times (on the order of 
days), although other distortions of the pattern may then 
appear. These are partially driven by a small drift of the 
pattern, which has a velocity on the order of one 
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FIG. 4. (a) Central section (43 x 53 cm) of a pearling state 
obtained by subjecting a pattern with wavelength 4.2 cm to a 
driving of a — 2.8 cm and / = ! Hz. (b) Corresponding power 
spectrum (range 188 x 94 cm~^) (c) Time series of PP when 
/ is changed from 2/3 Hz to /' — IHz and back again, (d) 
The order parameter APP as a function of the frequency /'. 
For details see text. 

ripple wavelength per 10^ oscillations The deforma- 
tions due to the drift are partly eliminated by using the 
slanted boundaries, allowing ripples to "drift out" of the 
system (leading to a small loss of sand at the edge of the 
plate). To obtain completely regular initial conditions, 
we have adopted the following procedure: First a flat bed 
is obtained by fluidizing the sand during a short period 
of strong oscillations of the plate. Then a regular pattern 
is imprinted into the sand by pressing down a frame with 
parallel equally spaced metal ridges. Small irregularities 
are then eliminated by a few 10) oscillations of the 
plate. 

By making large changes in a or / it is relatively simple 
to get regular secondary instability patterns as shown in 
Fig. ^ and |[ A study of the precise nature of these 
transitions involves runs performed at parameters close 
to the instability boundaries, and here the time scale for 
the development of instabilities diverges. As a result, the 
"ideal" dynamics may get hidden behind experimental 
artifacts. By introducing appropriate order parameters 
we will below precisely characterize both the pearling and 
the bulging instability. 

Pearling - When the driving frequency is increased 
beyond a certain critical value, we find a secondary 
"pearling" instability (Fig. Here, the crests of the ini- 
tial ripples remain essentially undisturbed, but in their 
troughs new small ripples (pearls) emerge, and the result- 
ing pearling pattern is periodic with the pearls aligned on 
lines inclined by approximately 45 degrees. The pearling 
transition is a supercritical secondary instability. When 
pearls are formed, their strength quickly saturates at 
some well defined value, and they disappear when the 
frequency is lowered as shown in Fig. He: Starting from 
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FIG. 5. (a) Central part (53 x 53 cm) of a bulging pat- 
tern, 525 T after the driving amplitude is changed to 4 cm 
(/ = 0.41). (b) The growth of deviations of a straight rip- 
ple pattern, quantified by the variance V for / = 0.55 Hz 
and a range of amplitudes. The thin curves with symbols 
correspond to amplitudes 2.4, 2.8, 3.0, 3.2 and 3.4 cm (below 
the bulging instability), while the thick curves correspond to 
3.6, 3.8, 4.0, 4.2, 4.4, 4.6, 4.8 and 5.0 cm (above the bulging in- 
stability) . (c) The maximum slope of V (for times up to SOOT) 
as a function of o shows a sharp transition for a ~ 3.55 cm. 

an imprinted perfect ripple pattern with wavelength 4.2 
cm (corresponding to an amplitude a — 2.8 cm), the sys- 
tem was driven at a low frequency of 0.67 Hz. After 1000 
oscillations, the frequency was suddenly increased to /', 
kept there for 500 oscillations and then again lowered to 
the initial frequency. 

To characterize the strength of the pearls, we have 
measured PP, the total intensity in the primary satel- 
lite peaks of the power spectrum (Fig. ^). Due to finite 
size effects and noise, PP is not zero for perfect patterns, 
and as order parameter we therefore use APP, the differ- 
ence between mean values of PP during period 1000-1500 
and period 0-1000. A plot of APP as a function of the 
quenching frequency /' shows a a well-defined transition 
point, above which the order parameter increases contin- 
uously (see Fig. IJd). A further increase of the frequency 
leads to more erratic states and finally to fiuidization 
where the ripple patterns are washed away. 

Bulging: ~ When the driving amplitude a is increased 
sufficiently, the regular ripple patterns become unstable 
to two dimensional modulations. While we cannot com- 
pletely rule out that this is a long wavelength instability 
0, our data strongly suggests that the wavelength of 
this modulation in the longitudinal direction is locked 
on four times the wavelength of the underlying pattern, 
with a similar wavelength in the transversal direction. In 
contrast to the pearling instability, this instability is sub- 
critical. The bulging deformations grow until neighbor- 
ing ripples become so close that they form defects, which 
climb and glide rapidly through the system, finally lead- 
ing to a regular pattern with a larger wavelength. For 
sequence of pictures showing the development of the in- 
stability see 

Close to the instability boundary it becomes very dif- 
ficult to distinguish slow drift from slow development of 
the instability, and we have therefore developed the fol- 
lowing sensitive measure for the onset of the instability. 
Starting with a perfect pattern of wavelength 4.2 cm, we 
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suddenly shift the amplitude. To characterize the time 
evolution of the pattern, we have extracted the local val- 
ues of the "longitudinal" ripple length Ai taken over the 
whole two-dimensional image. The variance V, defined 
as Il^i{Xi — X)^/{N — 1), is then a simple measure for 
the amount of deformations in the pattern. The evolu- 
tion of V is shown in Fig. for a variety of values of the 
amplitude a. 

Even below the formation of bulges, V grows slowly 
due to slow large scale deformations of the pattern, al- 
though the growth rate is essentially independent of 
the driving amplitude a. For a above some critical 
value, V displays a clear, surprisingly linear, increase 
with time. In Fig. ^ we have plotted the maximum 
of {V{tf + 100) - V{tf))/ 100 for t up to period 200. 
This quantity, which measures the maximum slope of V, 
clearly identifies the location of the secondary instability 
at a = 3.55(5) cm. In this way we can distinguish be- 
tween slow expansion or contraction of the pattern and 
the bulging instability. 

Doubling - When the amplitude a is decreased suffi- 
ciently, a sub-critical doubling instability occurs. The ini- 
tial phase of the development of this instability suggests 
that it can be captured in a one-dimensional framework. 
When the driving amplitude gets sufficiently small, the 
separation vortices that drive the sand transport do no 
longer reach over the through between ripples. This leads 
to the formation of bumps in the ripple troughs, which 
in turn grow out to form new ripples. Behavior similar 
to this has been seen in numerical studies of ID vortex 
ripple patterns P,p^ . For a picture of this transition see 

Discussion and outlook - We have characterized some 
of the pattern forming properties of fully developed vor- 
tex ripple patterns. By observing the number of moving 
grains on the ripples we have showed that the maximum 
shields parameter on the flat bed, calculated using the 
turbulent expression, is a relevant control parameter. We 
have shown that regular vortex ripple patterns are stable 
to small changes in the driving amplitude and frequency. 
When these changes become too large, however, the vor- 
tex ripples show a rich variety of secondary instabilities; 
pearling, bulging and doubling. 

Vortex ripples pose many theoretical challenges. The 
existence of a stable band and the doubling transition 
have been found in a simple model of ripple patterns 
The bulging and pearling are, to the best of our 
knowledge, not present in any simple theoretical models, 
such as Swift-Hohenberg type models incorporating local 
mass conservation, left right symmetry and a finite wave- 
length instability (which leads to spatial derivatives of 
order six). We believe that the origin of the instabilities 
is basically hydrodynamical and related to the dynam- 
ics of the separation zones. In the doubling transition 
it is thus clearly seen in the 1-d experiments that the 
new ripples originate approximately at the reconnection 



point for the separation vortex. The bulging and pearling 
transitions are genuinely 2-d and thus more complicated, 
but we speculate that the bulging transition is basically 
a Rayleigh-Plateau "sausage" instability of the almost 
cylindrical separation vortex, whereas the pearling insta- 
bility might be related to the centrifugal instability of the 
cylinder, giving rise to transverse Taylor vortices. Obvi- 
ously these ideas need considerable elaboration in view 
of the strong time dependent shear experienced by the 
separation vortex. 

It is a pleasure to acknowledge discussions with M.-L. 
Chabanol, J. Krug, A. Stegner and E. Wesfreid. 
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